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TECHNICAL NOTE 3768 


ON POSSIBLE SIMILAEIITY SOLUTIONS FOR THREE-DIMENSIONAL 
INCOMPRESSIBLE LAMINAR BOUNDARY LAYERS 
I - SIMILARITY WITH RESPECT TO STATIC®ARY RECTANGULAR COORDINAEPES 
By Arthur G. Hansen and Howard Z. Herzig 


SUMMARY 

Solutions of mainstream flow patterns for all possible incompres- 
sible lamlnar-bounda3:y-layer flows having classical similarity with 
respect to rectangular coo 2 rdinate systems are derived. These solutions, 
which apply to a wide range of flows, are summarized in table form. 


INTRODUCTION 

Turbomachine boundary layers are quite generally turbulent except 
perhaps in local regions. As illustrated by the secondary-flow behavior 
described in references 1 to 3, these boundary-layer flows are three- 
dimensional as well. Mathematically, the nonlinear partial differential 
equations that describe this turbulent three-dimensional boundary-layer 
flow are all but intractable. Nevertheless, it is important, for turbo- 
machine design procedures, to obtain theoretical solutions of the 
boundary-layer equations. In one much-used approach to this problem, by 
assuming two-dimensional laminar-boundary-layer flows, the equations are 
greatly simplified and solutions can be obtained. Actually, many useful 
applications of this two-dimensional laminar-boundary-layer flow theory 
have been made for estimating losses in turbomachines for design purposes. 
Often, however, this simplification is not physically acceptable. For 
example, when, in order to develop increased power output from compact 
engines, the flow velocities and mass -flow rates are Increased, the 
boundary-layer secondary-flow effects likewise are increased. The 
boundary-layer accumulations and blade -end-region losses associated with 
the secondary flows soon constitute a substantial portion of the turbo- 
machine losses. Unfortunately, these secondary-flow phenomena cannot be 
explained by extensions of the two-dimensional boundary-layer theory. 

In great measure, future progress in turbomachine design depends upon 
understanding and accounting for these secondary flows. As a consequence, 
the importance of understanding the three-dimensional turbulent-boundary- 
layer behavior has led to many experimental and theoretical investiga- 
tions in recent years. 
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Although siraplifying the hoimdary-layer equations hy two- 
dimensionalizing them appears untenable for high-output engines, one 
most important result of the experimental investigations of references 
1 to 3 is that laminar-boundary-3^gr behavior can usefully provide 
qualitative information concerning the turbulent -boundary-layer behavior. 
From a theoretical viewpoint, the importance of this observation can 
hardly be exaggerated. For the three-dimensional laminar -boundary-layer 
flows, the nonlinear partial differential equations are merely formidable; 
for the turbulent case the existence of turbulent fluctuation motions 
makes the determination of solutions much more complicated. The study 
of three-dimensional laminar-boundary-layer flows may, therefore, yield 
important information for the turbomachine designer. In addition, it is 
well to recognize that the extemial-flow wing-boundary-layer problem for 
high-altitude flight is also essentially a three-dimensional laminar- 
boundary-layer (secondary-flow) problem. (The existence of a limiting 
cross -channel-flow streamline depicted in references 1 and 2 eliminates 
consideration of a nonviscous-flow theory, at least in regions where it 
is reasonable to expect thin boundary layers . ) 

Consequently, theoretical investigations of three-dimensional 
laminar-boundary-layer flows over a surface having a leading edge have 
been made (refs. 4 to 8, e.g.). Two types of solutions are obtained: 
approximate solutions (perturbation method) for general mainstream flows 
confined to regions of small turning (refs. 4 and 7) and exact similarity 
solutions (for boundary-layer flows having affine velocity profiles) for 
restricted types of mainstream flows but with no restriction on turning 
(refs. 5, 6, and 8). The most general of these similarity solutions 
(ref. 8) is applicable for quite arbitrary main-flow streamline paths, 
provided the axial velocity component is constant (that is, constant 
static pressure axially) and the main-flow streamlines are all translates 
(no variation of streamline shapes in the tangential direction) . 

The present Investigation seeks to determine what further extensions 
of the similarity solutions are possible. Hew similarity parameters Ti 
are defined corresponding to mainstream flows not previously considered 
that enable the boundary-layer equations to be reduced to ordinary dif- 
ferential equations in terms of functions of these parameters. (Symibols 
are defined in appendix A.) These t)*s are found as a result of a 
systematic analysis using group theoretic approaches ("search for sym- 
metric solutions" and "inverse" method, ref. 9) described in appendix B 
of this report. The evaluation of the functions of tj and the deter- 
mination of the corresponding three-dimensional boundary-layer flow 
paths remain as problems of numerical analysis which are not considered 
here. Thus, while the integrations of the equations are not carried out, 
the similarity solution may be considered essentially complete for 
present purposes upon reduction of the boundary-layer equations to ordi- 
nary differential equations. 
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Solutions obtained Include cases of accelerating or decelerating 
main flows for quite general streamline paths, not confined to flow 
fields of streamline translates. These solutions are summarized in 
table I, which includes the mainstream velocity conrponents, the gen- 
eralized similarity parameter, and the final set of ordinary differ- 
ential equations for each family of mainstream flows. 

Similarity solutions obtained in previous investigations (refs. 4 
to 8 and 10 to 15) are likewise cited in the table to show their rela- 
tion to the solutions obtained in the present investigation. 


SBULAEITY SOLTJTIOIIS FOR GENERALIZED SDCLARITy PARAMETER ^ 

Similarity with Respect to Rectangular Coordinates 

The equations describing the steady incompressible thin-laminar- 
boundary-layer flow over a flat or sli^tly curved surface with coordi- 
nate axis oriented as shown in this sketch 



may be written 
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The houticLary conditions are 

u = V = w = 0 for y = 0 (at surface) 
lim u = U lim v = W 

y-> oo y 


(la) 


(lb) 

(l c) 


It is known from the results of references 4 to 8 that solutions involv- 
ing similarity parameters tj have "been bbtained. for the incompressible 
three-dimensional l ami nfl.r-'hf niTidp.-r y-l wy&r development on a nearly flat 
surface with a leading edge for certain restricted main flows. The pur- 
pose of the present investigation was to try to find solutions for much 
more general flows than had heen obtained heretofore by use of more 
generalized similarity parameters tj. As in appendix B and references 
4 to 8, the dimensionless velocity components u/U and w/w are defined 
in terms of functions F(t]) and 6(il)^ respectively, where tj ss T|(x,y,z). 
The continuity equation (Ic) may then be integrated, solving for v; and 
equations (la) and (lb) may be transformed (by substitution) to the new 
coordinates, x, z, and tj. The approach used herein is to determine the 
conditions und^ which the transformed equations (la) and (lb) become 
ordinary differential equations in F(t]) and G(tj), which can then be 
integrated numerically. For the purposes of this report, the similarity 
solutions will be considered completed when the defining equations have 
been reduced to ordinary differential equations. The actual numerical 
integrations will not be performed here. 


Generalized Similarity Parameter tj 

When the process outlined above is attempted, several practicsJ. 
restrictions appear at once. Unless the similarity parameter Tj(x,y, z) 
has y separable, that is. 


n = f(y)g(x,z) 

the transformed equations are inordinately complicated. Furthermore, in 
order to be able to integrate the continuity equation (ic) readily, tj 
must in addition be linear in y. Thus, tj is chosen as 
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( 3 ) 

and the integration can then he carried through. This definition of tj 
is not as restrictive as it first appears, because any similarity param- 
eter Tj* =Tj*(x,y,z) can he used for which a functional relation can he 
found such that 


fCl*) “ I =■ s(x,z) 

5ffiie functional relations of F and G to t]* are determined hy 

F = FCn) = F|f(Tl*3 
G = G(ii) = G|f(il*)] 

For an example, one could choose 

T]* _ 

(suggested in appendix B hy the similitude transformation) . The 
transformation 


= t, = y(x^z^-r)l/2r 

determines the corresponding parameter tj linear in y for use in the 
integrations . 

Thus any solutions obtained in the present investigation are clas- 
sical similarity solutions in the sense that, like all previous solutions 
involving a similarity parameter, the parameter is linear in y. 


Functions F(t]) and G(tj) and Boundary Conditions 

For flows over a nearly flat sxirface with thin boundary layers such 
as are being considered here, U and W are to he functions of x and 
z only, U = U(x,z), W = W(x,z). 

As indicated earlier, it is assumed that the velocity components 
u and w are expressible as 

u = U(x,z) F’(tj) (4a) 


and 
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w = ¥(x,z) G* (n) (4b) 

The stipulation that u and Mr vanish on the hounding surface and ap- 
proach the mainstream velocity at great distances from the boundary re- 
quires that 


F*(0) = G‘(0) = 0 

Um F* (il) = 1 
y-^“ 

11m G‘ (ij) = 1 

y-*.OB 

The form of the component v can now he determined from the con- 
tinuity equation. Substitution of equations (4a) and (4h) Into equation 
(ic) yields 

gf. +DF” ^ + + = 0 (5) 

Equation (5) can now he solved for v using equation (2), and the re- 
sulting solution Is given by 

(4c) 

For V to v ani sh on the hounding surface. It can he required without 
loss of generality that F(o) = G(0) = 0 and that f(x,z) s o. (Appen- 
dix C contains a discussion of the necessary and sufficient boundary 
conditions on F and G. ) 


Ordinary Differential Equations 

Substituting equations (2) and (4) Into equations (la) and (lb) 
followed by simplification then gives the following two dl^erentlal 
equations In f('h) and G('q): 

^ (p«2 _ _ 1) + w (G'F* - l) - || (GP") + 
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|| (G*2 - GG” - 1) + U <F‘G* - 1) - ^ (FG") + 

" (T) + " (^) - = 0 

As previously stated, the final goal of the similarity transformations 
is to convert the partial differential equations (l) into ordinary dif- 
ferential equations in F, G, and their derivatives only. Whether or 
not equations (6) and (7) become ordinary differential equations, how- 
ever, will ultimately depend upon the nature of U, W, and g2, as these 
quantities appear in the coefficients of all terms involving F, G, and 
their derivatives. The necessary requirement for retaining only terms 
in F, G, and their derivatives in equations (6) and (7) is that the non- 
constant coefficients of the various terms he proportional. If such 
proportionality exists, the coefficients can he divided out of the equa- 
tions, and a system of ordinary differential equations in tj will re- 
main. The condition for obtaining ordinary differential equations would 
then he equivalent to the following STOtem of partisG. differential equa- 
tions for the functions U, ¥, and g2; 


au 




a In U 


dW 


= a^lJ 


d m g^ 


dx 


_ d in g^ 2 TT o in 

= “5i 


( 8 ) 


where a^ 0 for i = 1, . . . , 6. The solutions of equation (S) 

describe all possible mainstream flows of the type U = U(x,z), W = W(x,z) 
for which solutions of the corresponding three-dimensional boundary-l^er 
equations could be obtained by the usual (i.e., linear in y) type of 
similarity transformation in rectangular coordinates. 


In the following section, certain consequences of requiring propor- 
tionality between the coefficients in equations (6) and (7) are investi- 
gated. Whenever such proportionality between two coefficients is assumed, 
it is referred to as an "ordinary-differential-equation condition" and is 
abbreviated "o.d.e. condition." (The o.d.e. conditions for the incompres- 
sible two-dimensional boundary-layer flows are discussed in refs. 15 to 
17. The o.d.e. conditions for the compressible two-dimensional boundary- 
layer flows are discussed in refs. 18 to 20 where systems of equations 
in two similarity parameters are obtained which closely parallel eqs. 

(6) and (7) of the present investigation.) 

Solutions for U, W, and g2 depend upon the o.d.e. conditions and 
the system of dlfferentleil equations in these quantities Trtiich results. 
However, the o.d.e. conditions in turn may depend upon certain basic 
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assumptions regarding the nature of U, W, and g2. For example, if it 

is assumed initially that U = U(x), the coefficient W — — will not 

appear in an o.d.e. condition. In general, it will he necessary to make 
such initial assumgptions about the form of at least one of the unknown 
functions before a unique set of b.d.e. conditions can be determined. 

The calculation of all functions can be carried out in a straight- 
f orward manner by making an initial assumption on the form of U or "W, 
by setting up the o.d.e. conditions, and by solving the corresponding 
differential equations. The cases where it is assumed that U = U(x,z) 
or W = W(x,z) are more complicated than cases resulting from other pos- 
sible assunqDtlons on the nature of II or W because a larger number of 
o.d.e. conditions result from equation (8) (if for no other reason). 

As an example of how calculation of the functions is actually carried 
out, the most complicated case, bhe case of initially assuming U = U(x,z), 
W = W(x,z), and g 0, is presented in appendix D. 


RESULTS AMD DISCUSSION 
Mainstream Configurations 

As a result of the analysis above, it can be shown that four families 
of solutions to equation (8) are obtained corresponding to 


( 9 ) 


u 

= me“^; 


w 

= ae^ 

u 

= mx?; 


w 

= pxF 

u 

= me“^z' 

r-1. 

N 

= ae“^! 

u 

= mx*^z^' 

-1. 

W 

= px"^"- 

to 

^U/5x 

in all 

cases) . 


the similarity between these resiilts and the results obtained for the 
two-dimensional case in reference 15 (pp. 116-120). For this case, the 
forms for U(x) (W is identically zero) are shown to be 


U = cx“ 


and 


U = ceP* 
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For convenience, the four basic solutions along vith five other 
forms (which, it turns out, correspond to the four basic forms under 
rotation of the coordinate axes) are presented in table I. OJable I, 
which includes within its organization all the possible solutions for 
mainstream flows having classical similarity with respect to' rectangular 
coordinates, for U = U(x,z) and W = W(x, z), is discussed in the follow- 
ing paragraphs . ^ 

U and W . - U and W, the velocity component^ of 'the mainstream 
in the x and z directions, respectively, are listed first in the 
table. 

For all except case VI (obtained from case IV by rotation), the 
chart indicates U and W are representable by (a) powers of linear 
functions in x or z individually, (b) exponential functions in x 
or z alone, or (c) the products of such functions. In case VI, U and 
W take the fonns of powers of linear conibinations of x and z result- 
ing in mixed polynomials in x and z. It may be observed that great 
freedom of choice of U and W velocity components can be obtained by 
assigning different values to the available constants. 

Erojection of mainstream on surface . - The equation for the projec- 
tion of the mainstream on the surface may be obtained by integrating 
dz/dx = W/U. For flows nearly parallel to the surface, this projection 
approximates the actual main-flow streamline path. When the constants 
are chosen such that straight mainstream flows result (always true for 
cases III and VI ), there can be no boundary-layer crossflow relative to 
the mainstream. For such flows, F = G and equations (6) and (7) are 
identical. It may be observed that, except for solutions cited under 
case IV and the stagnation solution of reference 10, only solutions hav- 
ing straight streamlines (in the x-z plane) have been completed. 

Irrotationality . - This listing indicates choices of constants for 
which the mainstream component of vorticity normal to the STirfewie 
(Su/^z - 3w/Sx) is zero. For the flows considered here, only this com- 
ponent of vorticity can be much different from zero, in any event. 
Therefore, specifying dU/Sz - Sw/Sx = 0, as is done for the listing of 
irrotationality, actually serves to set the conditions for nearly irro- 
tational main flow. 

Irrotationality and two-dimensional continuity . - The flows repre- 
sented by cases III and VI are straight flows. Cfbviously these flows 

can be irrotational and have two-dimensional continuity ^ 

only for the case of constant velocity. 0?he combined conditions would 
require strai^t constant-speed flows for cases I, II, IV, V, VII, and 
VIII as well. Only case IX can satisfy both irrotationality and two- 
dimensional continuity for nonstraight flow paths, when n = 1, r = 1, 
and m = - p. The main-flow streamline paths here are hyperbolic. There 




10 


MCA TN 3768 


are, however, difficulties in inteipretation of the significance of case 
IX flows and. their associated, boundary layers which will he discussed 
later. 


Boundary Layer 

SiTiHiarlty parameter ii . - The so-called similarity parameter de- 
fined as T) = y g(x,z)/>\/^ provides information concerning the physical 
dimensions of the boundary layer. At the surface = 0, and as 

the boundary layer merges into the mainstream. However, some finite 
value of 'H = 'Ho serves as a practical outer limit to the boundary layer 
(refs. 4: to 8). In terms of ph^ical lengths, the boundary-layer thick- 
ness at any point is the hei^t y at irfiich the value of reaches 
approximately tIq. From y = outer edge of the 

boundary layer, it can be seen that the boundary-layer thickness is in- 
versely proportional to the magnitude of g at any point on a surface. 

Thus, in case IH, for example, for n and c positive integers, 
the boundary layer becomes progressively thinner with increasing x and 
z because g(x, z) becomes progressively larger. 

g(x,z) constant . - When g is constant, the boundary-layer thick- 
ness does not vary over the surface. This corresponds, for example, to 
the stagnation-point flows (refs. 11 and 15). 

Eeference 15 notes that in the two-dimensional stagnation-point 
flow with the constant boundary-layer thickness (g constant) an exact 
solution for the Havier-Stokes equations is obtained where U = ax and 
V = ay, for constant a. It is interesting to note that the three- 
dimensional stagnation-point flow of reference 10 (case IX) with constant 
g likewise is an exact similarity solution for the three-dimensional 
Havier-Stokes equations with U = ax, ¥ = cz, V = - y(a + c), g^ = b, and 
T) = y>\/b/v. For this case, however, the Havier-Stokes equations actually 
reduce to the boundary-layer equations without casting out terms because 
of physical considerations, such as very thin boundary layers, and so 
forth. 

Le a din g edge. - The problem of relating the mathematical solutions 
to physical reality is Intricately involved with the nature of and 
g(x,z). Usually the boundary layer must be considered as "beginning" 
somewhere in a real flow at a leading edge. At the leading edge the 
boundary layer haa zero thickness and develops under the main flow along 
the surface. A line on the surface along which t] takes on very large 
vailues (because g(x,z )-►“’) maybe considered such a leading edge. Along 
such a line the boundary layer has zero thickness being inversely pro- 
portional to g(x,z) . Because lim F’ = lim G’ = 1, then u = UF’ = U 
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and w = WG' = W and the houndary-layer velocity components merge smooth- 
ly into the corresponding mainstream components. An example occurs in 
case IV with n = 0 along the line x = 0 on the surface. ODiis, of 
course, is the case treated in reference 8 as well as elsewhere. 

0?here are situations where the mainstream velocity components are 
zero along a line or along an entire plane normal to the surface. It 
is readily seen that the houndary-layer velocity components are also 
zero matching the mainstream. If cases of g constant are not included 
for the moment, inspection of the chart reveals that either .g(x, z) like- 
wise becomes zero (case IV, n>l along the line x = O) and so the 
’'boundary layer" is infinitely thick along such a line or g(x,z) becomes 
very large and the boundary layer begins at the line (case VI, n = l/2 
along the line ax + cz = O). The situation where g(x,z) = 0 along a 
line is typical of this type of accelerated flows with a bound^y layer 
becoming progressively thinner along the surface. In such cases,, the 
sinrplest thing to do is to confine the discussion of the flows to such 
regions where the requirement of thin boundary layers. is satisfied. OSie 
correspondence in these canes between the mathematical solutions in the 
thin-boundary-layer region and physical reality can be properly deter- 
mined only by eaperlment. 5Phis, however, is equally true for flows with 
satisfactorily defined leading edges. 

When g is constant as in the stagnation-point flows discussed 
previously, again no leading edge providing a "beginning" place for the 
boundary layer can be defined. Likewise, the physical realism of a 
solution is questionable when the "leading edge" is at x or z » * or 
when the velocities at an otherwise well-defined leading edge (tj = ») 
take on infinite values. Case IV when n = 0, U = Uq., and as 

X-+-0 presents a typical case of flow over a surface with a leading 
edge. When r<0 however, in this case (refs. 7 and 15, e.g.) the 
tangential component W takes on infinite values along the leading 
edge. 


Ordinary Differential Equations 

The actual numerical solutions of the ordinary differential equa- 
tions are beyond the scope of the present investigation. Numerous ex- 
amples have been calculated elsewhere for particular values of the free 
constants in the equations. Many of these calculations are Indicated in 
the listing "References and comments" associated with each cane. While 
the existence of a solution and its calculation must remain to be deter- 
mined individually for each case in general, certain remarks can be made 
here concerning particular situations. 

Separation of F and G . - Numerical solutions will be considerably 
easier in those cases for which only one of the parameters, F or G, and 
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its cLerivatives, appear in one of the two equations. The listing "Sepa- 
ration of F and G" enumerates conditions where this will occur. Under 
these conditions the ordinary diffepential equations reduce to the forms 
of cases I and V. 

In case I, equation (6) is a Faliner-Skan equation (ref. 13) as is 
equation (7) in case V. In all cases, for choices of the constants that 
effect separation of F and G in one of the equations (conditions 
which are most likely to he useful for purposes of computation) the nu- 
merical solutions for the equations in F alone or G slone are already 
known (refs. 14 and 15). Likewise the numerical solutions for equation 
(6) (cases IV and Vl) and equation (7) (case H) are known. 

Althoiigh it is not apparent in the table, if a = 0 in cases I, 

III, and VII, or p = 0 in the remaining cases, then w = 0 and equa- 
tion (lb) and, hence, equation s(7) disappear. Under these con(^tions, ex- 
cept in cases II and V (which are simply rotations of cases I and IV, re- 
spectively), equation (6) becomes the FaJkner-Skan two-dimensional flow 
equation for which the complete solution is known (ref. 13). Sitmilarly, 
whenever m = 0, u = 0, equations (la) and (6) disappear, and equa- 
tion (7) becomes the FaHmer-Skan equation except in cases I and IV. 

Boundary conditions not achievable . - Should the ordinary dif fer- 
ential equations reduce to the form =0 (or G"’ =0) in any case, 
the solution is 


F = kiTjS + kgTJ + kg 

However, the boundary condition lim F*(ii) =1 cannot be achieved. The 

Tj — >■ oo 

occurrence of such forms of the equations is listed as "Boundary condi- 
tions not achievable." The requirement that F’ and G* approach their 
limiting values asymptotically as t) ->oo restricts somewhat the choices 
of the constants as well (ref. 19). For example, when F* is to ap- 
proach its limit 1 from below, the relations must be such that the cur- 
vature of the F* curve, which depends upon F"* , must be negative. 

That is, from equation (6) 

j-nt = F"t (10) 

The practical analysis is somewhat simplified for conditions when separa- 
tion of the variables occurs. Then, F* is the Falkner-Skan function. 

F, F*, and F" are positive, and cheeks can be made readily during the 
numerical calculation procedures for the appropriateness of the sign of 

F..I . 
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Similar statements can l)e made for the restrictions on the constants 
required for correct valties of G”* and for conditions when the functions 
approach the limit 1 from above. 

Linearity in u or v . - Under conditions for which there is line- 
ar ity’~oF'eqiIatinnB~"(l5"'i^ u or w, that is, ha/hx e o or bw/hz h o, 
an extension of the solution heyond strict similarity of the velocity 
component with respect to the co 2 nresponding x or z coordinate can he 
made. Reference 8 under case VII makes use of such linearity with 
respect to w to obtain solutions for quite general streamline shapes. 

In order to simplify the computations, U = Uq (the cane when n = O) was 

chosen in reference 8, so that F becomes the already well-tabulated 
Blasius parameter F, (ref. 15) leaving only the parameter G to be 
computed. 


CONCLUDING REMARKS 

Solutions for mainstream U and N flow components for all pos- 
sible boundary-layer flows having classical similarity (having affine 
velocity profiles) with respect to rectangular-coordinate systems are 
obtained in an orderly fashion. QJhese solutions, which apply to a wide 
range of flows, are summarized in table I. Careful attention should be 
paid to certain practical considerations. 

The final solutions for the boundary-layer velocity components have 
not been carried out. The actual numerical solutions for the ordinary 
differential equations derived here are considered beyond the scope of 
this Investigation. They belong more properly to a program involving 
the use of high-speed conputing machinery. 

It is important to note that there is really no case described for 
such simple configurations as flows through a row of stator blades with 
appreciable tviming, irrotationality, and very near two-dimensionality 
of the mainstream. The cases for flows with a satisfactorily defined 
leading edge include but a portion of the totality of flows described. 
Nevertheless, physically, the beginnings of a boundary-layer development 
would reasonably seem to be most important. Except for that described 
in reference 8, there is little experimental verification available for 
any three-dimensioneil simillarity solutions and none for most solutions 
indicated in table I. 

On the positive side, there are obviously conditions for which the 
solutions presented enable reasonable approximations to be calculated 
for three-dimensloneJ. laminar -boundaiy-layer behavior (see ref. 8) . More 
Important, however, if no case can be found in table I to serve as a 
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reasonalDle approximation for a particular flow problem, then in all 
likelihood some approach other than the similarity solutions in rectan- 
gular coordinates must he sought. 


Lewis Fli^t Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, June 6, 1956 
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AEPEHDIX A 


N 

O 

O 

lO 


symbols' 

constants 

f arbitral:^ function 

function of similarity parameter, u sUP'Ctj) 
function of similarity parameter, v s WS'(t]) 
g,g(x,z) function of coordinates, x z 


h,k,m,n, 

U,V,W 

u,v,w 

C,v,^ 

x,y,z 

x,y,z 

a 

n 

^0 

Tl* 

V 


constants 

mainstream velocity conponents in x,y,z directions, 
respectively 

boundary- layer velocity components in x,y,z directions, 
re^ectlvely 

rectangular coordinates 


constant 

simil a rity ( space) variable, t] = g(x, z) 

value of T] 'Which defines practical outer limit of boundary 
layer 

similarity variable, tj* = Ti*(x,y,z) 
coefficient of kinematic viscosity 


Subscripts: 

0 constant 


1 = 1 , 2 , 

3, • • 


index numbers 


Superscripts: 

Primes denote differentiation 
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APPENDIX B 


SEAECH FOR SOLUTIONS 
Approach, to the Problem 

One puipose of the present investigation was to develop a logical 
systematic method for obtaining more genei^al solutions for three- 
dimensional Inconpressible thin-boimda^ry-layer flows having affine dimen- 
sionless velocity profiles (so-called similarity solutions) then hereto- 
fore available. The utility of In^ectional analysis as a method for 
obtai n i n g such solutions to the partial differentieil equations of fluid 
flow is discussed in chapter IV of reference 9. The basis of the method 
is the axiom (obtained from group theory) that, if the hypotheses of a 
theory are invariant under a group of transformations, then so are its 
conclusions. Si:q)pose a group of transformations can be found that leaves 
the form of a system of partial differential equations invariant. The 
method, called the "search for symmetric solutions", as applied to 
boundary-layer theory involves finding combinations of the reference 
coordinates that are likewise Invariant under the groi^. Of even greater 
significance to the present investigation is the concept of an "inverse 
method" (ref. 9, p. 140), which states that, given a partic ular solution 
with certain properties, it may be possible to find additional solutions 
by a priori postulation of those properties. The procedures for applying 
these two methods to boundary-layer problems are described in detail 
later. It appears advisable first to establish the reasonableness of 
anticipating that the application of these methods can lead to solutions 
of the boundary- layer equations. 

Reference 9 applies the "search for symmetric solutions" to the two- 
dimensional laminar -boundary- layer equations in rectangular coordinates 
and is able to derive a systematic, logical development of the Blasius 
solution in terms of a similarity parameter t] = y/V^* Furthermore, it 
is known from the results of references 4 to 8 that solutions involving 
such similarity parameters have been obtained for the incon^ressible 
three-dimensional 1 aminar-bovmdary-layer development on a nearly flat 
surface with a leading edge for certain restricted main flows. Thus it 
seems reasonable to attenpt to apply the method of "search for symmetric 
solutions" and the "inverse" method in order to find three-dimensional 
boundary- layer solutions for more general main flows. 

G?he procedure for this investigation was to use the "search for 
symmetric solutions" in an effort (l) to find some solutions for the 
incompressible three-dimensional laminar-boundary- layer flow over a sur- 
face in a logical, systematic manner, (2) to study the properties of these 
solutions, and (3), if possible, to invoke the inverse method to find 
other related solutions. In this fashion, it weis Intended to try to 
find solutions for much more general flows than had been obtained 
heretofore . 
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Following reference 9, under transformations of similitude (changes 
of scale). 


. X -»■ a^x, y -*■ ttgy. 

Z ->• ttgZ '> 

u -♦ a^S, V -*■ OgV, 

W -»■ OqV \ 

U -*• OyU, V Oq V, 

W dgW J 

the system of equations (l) is invariant when. 

for example 

_ _ — —2 _ ~2 



_ -1. _ 1 


► 

“S “ °'8 “ ~ a- 


II 

II 

II 

s 

II 

1 

J 



(Bl) 


(B2) 


where the last constants were chosen equal to 1 for convenience. The 
invariance of equations (l) under this transformation may he verified hy 
direct substitution of equations (Bl) and (B2) into equations {l) . 


As mentioned previously, the first step in the search for symmetric 
solutions is finding combinations i] of the reference coordinates that 
are likewise invariant under the transformation of similitude using equa- 
tions (Bl) and (B2) . Several such combinations suggest themselves at 

once: tj ~ y/ ^/x + a, y//x*^“°-“^, y^/ + z^, and others. 

Substitution from equations (Bl) and (B2) verifies the Invariance. The 
solutions are to be obtained by expressing dimensionless u and w as 
functions of q under conditions that will enable these functions to be 
evaluated. By conparison with reference 9, it can be shown that solutions 
which may be obtained in this fashion will be similarity- type solutions 
with more generalized ^ace variables tj than considered heretofore 
(refs. 5, 6, and 8) . 


Application of Method 

Specifying the dimensionless velocity components u/U = F'(ri), 
w/W *= G'(t]) as in equations (4a) and (4b) of the text and defining 

T] = y/yy/x + z, for an example, it is found that the transformed boundary- 
layer equations reduce to ordinary differential equations in P(ti) , 

G(ti), and their derivatives for 


U » m 
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W = p 

V = 0 

This is the case of a family of straight main flows (actually a 
rotated Blaslus flow) with streamline projections in the x-z plane de- 
scribed by z = (p/m) X + Zq. 

Similarly, as suggested earlier, when 

T] a y/^2c“z~^“^ 

the transformed equations reduce to ordinary differential equations for 

U = m(x)^‘*'^(z)““^“^ m / 0 

W =p(x)^(z)~^ P / 0 

Here is the family of mainstream flows with projected streamlines in the 
x-z plane described by z°^ = ZqX®. 

Thus, preliminary application of the method of "search for symmetric 
solutions" using transformations of similitude has led directly to ob- 
taining two new similarity solutions. In order to invoke the inverse 
method- and perhaps extend these results, the solutions were studied. As 
a result it was found conveaaient to postulate a priori the property that 
the boundary-layer equations can be transformed to equations containing 
functions of a generalized similarity parameter t], and that the trans- 
formed equations can be reduced to ordinary differential equations by 
determining specific forms U(x,z), W(x,z), and ii. Four families of 
solutions are obtained upon application of the Inverse method in this in- 
vestigation and are listed as nine general cases in table I. 

The associated problem of actually determining the numerical solu- 
tions to these families of differential equations is twofold. There is 
first the problem of the existence of such solutions for any particular 
case. As discussed in the text, values can be diosen for certain free 
constants in the equations such that the required boundary conditions 
cannot be achieved. Second, there is the fact that, for even moderate 
ranges of the constants that enter the equations as parameters, literally 
thousands of essayable Integrations are represented in the table. The 
scope anfl purpose of this investigation is to present the flow equations 
in such form that their applicability to the boundary- layer development 
under some particular mainstream configuration will be readily apparent. 
The problems of numerical integration for that particular case can then 
be undertaken more appropriately at that time (as in refs. 4 to 8 and 
10 to 15) . 
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APPEHDIX C 
BOUICDARr OOHDITIONS 

The houncLary- layer, equations (l) have "been transformed using the 
definitions 


u = UP'(ti) 

(4a) 

*! 

Ill 

B 

(4b) 


and the solution for v is given hy equation (4c) . 


The mainstream flows have then been determined for which the trans- 
formed boundary- layer equations reduce to ordinary differential equations, 
involving F(ii), G(ti), and their derivatives. The boundary conditions on 
these functions of t) must be chosen to correspond to the boundary con- 
ditions on u, V, and w in' the original partial differential equations 
(l) . Corresponding to the four boundary conditions, 

u = w = 0 for y = 0 


lim u = U 
y-4>» 


lim w = W 

W‘-»- • 


the following four boundary conditions on P and G result from equa- 
tions (4a) and (4b) 

P«(0) = G*(0) = 0 
lim F‘ = lim G* « 1 

T)-> OB TJ-*’ OB 

The remaining boundary condition on equations (l), v = 0 for y = 0, 
requires that (see eq. (4c)) 

(i - 1 - (n - 1 

Substituting equations (4) into equations (la) and (lb) with f(x,z) 
defined as in (Cl) yields 
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^ ^ (G’F* - 1) - H P"[G - G(0)] + 

U El tF - p(0)] + ^ il [G - G(0) ] - g2p...^ 0 

ox 2 oz 2 

(C2) 

^ |g‘^ - G" [G - G(0) 1 - :^+ ^ ^ (I^'G’ - 1) - ^ G" p - P(0)] + 

^ £(j _ Q(0)] + U - — [F - F(0) 1 - g^G"’ = 0 

(C3) 


If the substitution is now made, 

F = F + F(0) 

G = G + G(0) 

!Phe above equations reduce to the form^^iven in equations (6) and (?) 
with F and G replaced by ? and G, reflectively. Provided that the 
coefficients of (C2) and (C3) have been made proportional, the solutions 
for equations (C2) and (C3) are uniquely determined by solutions of equa- 
tions (6) and (?) in P and S withjjhe boundary conditions 
fCo) a S(0) a F'(0) = G’(0) = Oj lim F* = lim. 5’ = 1. Furthermore, 

eo oo 

U = UF' a UF' (C4) 

w=WG*aWG» (C5) 

F(o)l + H fG - G(o)3} + {nF* - Lf - F(0)]j. + 

- [G - G(0)]} (C6) 

F + || G^ + ^ (t]F« - F) + W (tjG* - G) (C7) 

Here the solution for u, v, and w is^uniquely detennined by the solu- 
tion for F and 5 . Finally, F and G can be identified with F and 
G in equation (4c) and the provision f(x,z) a 0 in (4c) and 
F(0) a G(0) a 0. 
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AEPEKDIX D 


EXAMPLE OF DERIVATION OF EXPLICIT FORMS OF U, W, AND g2 

The cases of obtaining the ordinary differential equations for flows 
ass umin g U = U(x,z) and W = W(x,z) are presented here because they 
Include all the variety of situations and conplicatlons that arise with 
other possible assun^tions on the nature of U wnfl This assuB^tlon 
is equivalent to the hypothesis that the first partial derivations of U 
and W do not vanish Identically. For convenience in the analysis the 
seven different coefficients (see eq. (8)) are listed in the order of 
their appearance in equations (6) (7): 

® au/Sx 
(D hM/hz 

© 

As Su/3x and Sw/dz have been assmned to be nonzero, an o.d.e. condi- 
tion for @ and © will be 


Hence, 


Similarly, from (|) and © 


au „ d In W 

1 SE - ° -SE” 


b. 

W = fi(z) U 

(Dl) 

W = fgCx) 

(D2) 


Now, if bi = b2, it follows from equations (Dl) and (D2) that fi(z) 
and f2(x) constants. Hence, 


W = bgU 1 


(D3) 
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If "b 2 ^ then from equations (di) and (D2) and @ @ 



U = f3(x)f4(z) 

(D4) 


W = h^f^(x)fjz) 

(D5) 

hi ^ h 2 " - First consider the case h^^ / ^ 2 - 

Then hy @ and (§) 


g^ * hgf^(x)fj^(z) 

(D6) 

As Bw/Bx = h.f 3 (x)f.(z) 
then 

from (D5), it follows that, since aw/ax 0, 


a In g2 _ j_ 

~SJ — 

(D7) 

Similarly, Su/Sz 0 in^lles that 


a In g2/az = f^(z)/f^(z) 0 

(D8) 

Hence, from (X) and 

!i ^ £i 

f ** 6 fi 

^3 ^3 

(D9) 

or 

-L 

■^6 

= V3 

(DIO) 

If hg = 1, 

h„x 

(Dll) 

If hg 1, 

h,, 

^3 == (^9^ + ^10^ 

(D12) 


The constant represents only a displacement of the coordinate axis 

a.nfl may he taken equal to zero without loss of generality. 
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Similarly, from @ and © , and equations (D5) and (D6) 


^4 - ^ 12 ® 


\3" 


(D13) 


or 


^4 = ^ 14 == 


15 


(D14) 


From the various combinations of f3 and f^., the following forms 
for II and W result: 


U = me^e°="j 
U = me°^^-lj 

U a mx“e“5 

TT — ^"1 

U = mx z 5 


.. nx cz 
W = ae e 


W a ae^^ 

TT n-1 cz 
W = px e 




W = px'^'^z^ ^ 


(D15) 


hi a h2* - Now consider the remaining case h]^ a h2» If "bl = 0^ It 
follows from equation (D3) that SW/Sx a Sw/Sz a 0, which contradicts the 
basic hypothesis. If h^L ® ^ “ ^3X1. Then from @ and (f). 


If fc 0, then from @ and 



(D16) 

U a fg(c^X + z) 

(D17) 

w a "bgfgCc^x + z) 

(D18) 

independent variable. 

f ^ a 0 inqplles that 

fg a CgCc^X + z) 

(D19) 

U a CgCCjX + z) 

(D20) 

W a bjCgCc^x + z) 
and © 

(D21) 


0 (D22) 
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Similarly, from (|) and 


« V -pii i n 

35 " “ — 2 = f 0 

oz 


Therefore, from @ and @, 


°1^5 = °5°1 f . 

5 


^5^5 


and therefore. 


If c. = 1, 


^5 - 


fg = c^e 


Cg(Cj^X + z) 


If C 5 


(D23) 


(D24) 

(D25) 

(D26) 

(D26) 


- t=g( 


=2^ + 


“) + “loJ 


'11 


(D27) 


As "before, c^q taken equal to zero without loss of generality. 

.The same result is obtained by use of (D anfl . 


In summary, for b]^ = 1 the following forms for U 


and W result: 


TT cz 

U *= me e } 

U = m(ax + z)’^; 


TT nx cz 
W = ae e 

W = p(ax + z) 



(D28) 


The final case to be considered is b^^ 1. For b]^ 1, cannot be 

identically constant . If = k, a constant, then from Qi w-nd (§) U 
would have to be of the form U = d^x + h-^Cz). Then from @, wnd 

equation (D3), b 3 [dix + h]^(z) ]^”Ih^(z) = dg. But this requires that 
either b^ = 1 or h4(z) =0* Neither case is admissible. Furthermore 
neither og/Bx nor dg/dz can be identically zero if the hypothesis 
that the first partieil derivatives of W and U do not vanish is not 
to be violated. If for example g2 ss g^(x), then from @ and ©it 
would follow that U = U(x) , and hence SU/dz = 0 . 
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i 

© 

@ the following form for U is determined: 



u = (g2) ^2^2;) 

(D29) 

From @ and © , 

W = (g2) ^g(x) 

(D30) 

However, 

W = hgU^^ 

(D3) 

Therefore, 







iigCx) 


Now if ^4 / ^1^3^ follows at once that U wnfl w can he e^garessed 
as a product of a function of x arid a function of z. Hence, equation 
(D4) would hold, and the analysis would proceed accordingly as before. 


If d 4 = h 2 _d 3 , then hgCz) and h 3 (x) must he constants and ®, 
(§), and equation (D29) give 

2 3 . 

dgCg^) 2 = SgV^x (D31) 

From ©, and equation (D30), 



(D32) 


If dg a 2 or d^ = 2, then g2 = g(x) or g^ = g(z); and the hypothesis 
is violated. If dg = 1, then from equation (D3l), 


a In g^ 

5x 




(D33) 


= h4(z)e^^ 


(D34) 


Similarly, if d^. *= 1, it follows from equation (D32) that 

g2 = h5(x)e^^ 


(D35) 
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If d .3 1,2, integration of (D3l) yields 

1 

^ 1 

+ Ii 6 (z)(d 3 - 1)] (D36) 

If d^ »= 1,2, integration of (D32) yields 

1 

g2 = [{dgz + h 7 (x))(d 4 -1)] (D37) 

Eq,uations (D34), (D35), (D36), (D29), and (S30) require that U and W 
he one of the forms already given in (D15) or (D28) . 

QThese forms thus obtained constitute a-H the possible solutions 
etrisj.ng from the initial assunptlons . 
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TABLE I, ~ SDaLABITT BOLOSima IN STMIONART HECEAWDIJR OOCmONiZES 



Oaae I 

Oaae II 

Oaae III 

n 

ue»“ 

ina>u 

Bw^e°» 

V 

aa** 

pa“ 

ee"*e«« 


^ -AtAV 2 _ r w , 

_Ae««\i/a r OT *iV2 _(^V2 

< 

l>oiWe<»«\Va _/bUNV2 



V ) ° " y{ai/ 

Bquatloa (6)' 

im[{P')2 - ^ - ij- bF"' . 0 

-op np(0'F' - 1) - UFin x 0 

1 

J 

(i»fan)[(Fi)2 - I|^ - ij- bF'" - 0 
1 

F - 0 

Equation (7) 

■r(F>a> - 1) - nn ^ - ba»< •• 0 

op[(Q')® - ^ " ° 

Projection of main 
atrean on turfaoa 

■^r-n^ ® '*• “0 

*'’pI§=oT-^"-”^’ + *o 

a a ^ *0 

Irrotatlonallty 

ar » 0 

Bin X 0 

on " an 

Linearity In u 

nn a 0 

Linear 

nn - 0 

Linearity In w 

Linear 

op X 0 

ao X 0 

Boundary ooodltlona 
not aohlevable [t1 

mn ■ 0 
DB> ■ nm ■ 0 

op X np X 0 
op X 0 

j Bn X -CO 

Separation of (6} 

F and a (7} 

Separated 
■1 > nr - 0 

op X up X 0 
Separated 

1 Soaa not apply 

Referenoea and 
oocmantB 

Qeneral aolutlon of eq. (6) in refa. IS 
and 14. 

If n > r, then F - G. 

Oaae I with rotation of coordinate axle. 
Eq. (7) la aolTed in refa. IS end 14. 

If 0 X n, F X a, snd eq. {6} la oon- 
pletely aolved tn refa . 13 and 14. 

Thla la a oaae of rotated two- 
dlaenalonal flow aolved tn refa. 
IS and 14 with b - an + ao. 
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TABLE I, - Concluded. SIMILARITY SOLUTIONS IN STATIONARY RBOTANGULAR COORDINATES 



Oasa 711 

Onao TOt 

Oaao IX 

u 

M»V-1 


■nV-l 

V 


p^i«“ 

paP-V 


AonylVA _ ^\l/4 _ (M yA 

A^odi/fe _ ^,\Vj _ 

^ \i/a Aw \i/ii 




Iqnablon {() 

«a|^{P')®-i5l-l]+a(r-l)(0iP'-l)-«{i+l) Bp-bPin-O 

■iIjtF'}‘-l]-a(nfl) 2f+oi»(aipi-^-:^-bF«n -0 

Ml|][F')*>3-«(Dfl) + 

p(r-l)(a<F'-]) -p(rbl) flp.-hF'« - 0 

Bq^iaUoa (7l 

ojf -^-W™ - 0 

op[(a')’ - ^ ^j4B(n-i) (r>ai-i) .a(iM-l) -bO"' > 0 

P» Ijo ' )®- l] -p(rt-l) + 

b(d-I) -aCttfl) -tadi* >0 

hvjartlon or atia 
■tnia on nrCiioa 

■ - - In « + ^) 

X - - In X + 31^,) 

«o*" ■ ^ 

InrataUonolltr 

n(r-l) ■ sn • 0 

BO a p(n-l) >■ 0 

n(p~l) " p(n-l) ■ a 

LUMorttr In u 



n ■ 0 

Unoanor In v 

a « Q 

«P H 0 

p - 0 

llilgfirptffi HilMnl 



n ■ p ■ 0 

a - p " 0 

a«parttlon or (6) 

F and 0 ( 7 ) 

a « 0 

CBS ■ D 

op M 0 
n « 0 


Rararonaoe and 
oeoanta 

If ■ - 0 , tbra MO. (In) nd (•) dloapiionp, ' 
p D Q, and M. 17) la ooaplntalp aolvod in rafa. 
la aid 14. , . 

ir a > 0 , Ww M>. (lb) and (7) dlaamar. a ■ 0, 
and oq. (6) la aolvad in rafo. IS and 14. 

Bm oaia 7Sa 

If p - 0, aqa. (l^) and (7) 
diaappaap^ Q n 0 . and m. (6) la 
solTid In rafoi IS and 14# 

If ■ ■ 0| ovo (iQ and (6) 
dlaappoorj F n 0 . and ta. (?) la 
lol'nil In nafa. 13 and 14. 

If p " n, F .. b, raf. 13. 

Ref, lOi ran.nabjpaa (nf. 10). 

Bof. Ill plana abasutlon pMM.nowi 
n > 0. n ■ If B - b. In M. W. 

Bar. 13| p. 71i plana ataBnatlcn flow: 
a • 1. a - B| p - 0. 

» 










































